Hyperdiffusion in non-linear, large and small-scale turbulent dynamos 
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The generation of large-scale magnetic fields is generically accompanied by the more rapid growth 
of small-scale fields. The growing Lorentz force due to these fields back reacts on the turbulence 
to saturate the mean-field and small-scale dynamos. For the mean-field dynamo, in a quasi-linear 
treatment of this saturation, it is generally thought that, while the alpha-effect gets renormalised and 
suppressed by non-linear effects, the turbulent diffusion is left unchanged. We show here that this is 
not true and the effect of the Lorentz forces, is also to generate additional non-linear hyperdiffusion 
of the mean field. A combination of such non-linear hyperdiffusion with diffusion at small scales, 
also arises in a similar treatment of small-scale dynamos, and is crucial to understand its saturation. 
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Introduction: Large-scale magnetic fields in astrophys- 
ical bodies are thought to be generated by dynamo ac- 
tion involving helical turbulence and rotational shear Q . 
For turbulent motions, with a large enough magnetic 
Reynolds number (R m henceforth), this is also accom- 
panied initially by the more rapid growth of small-scale 
fields, correlated on the tubulent eddy scales and smaller 
d H 0. An important problem is to understand how 
the Lorentz forces due to these fields back-react on the 
turbulence and hence lead to mean-field and small-scale 
dynamo saturation. 

Semi-analytic treatments of the back-reaction have 
typically used the quasi-linear approximation (see below) 
or closure schemes to derive corrections to the mean-field 
dynamo coefficients. It is then found that the a-effect 
gets "renormalised" by the addition of a term propor- 
tional to the current helicity of the small scale fields. 
But at the same time the mean field turbulent diffusion 
does not get affected, if one imposes the incompressibil- 
ity condition on the velocity field (including the compo- 
nent induced by the Lorentz force) 0, |f| . This is perhaps 
somewhat intriguing, as one would have naively expected 
the Lorentz forces to affect all the transport coefficients. 
Further, if the effective "turbulent diffusion" does not get 
modified at all due to non-linear effects, one also wonders 
how the non-helical, small-scale dynamo would saturate 
at all? We clarify these two issues here. 

Large-scale dynamo: The induction equation for the 
magnetic field is given by, 



<9B 

— = Vx(vxB-iiVxB) 
dt v ' ; 



(1) 



where B is the magnetic field, v the velocity of the fluid, 
and r\ the ohmic resistivity. In the kinematic limit, it is 
usual to take v = vo + vy, the sum of an externally pre- 
scribed large scale velocity field v and a random field vy. 
Also vj- is generally assumed to be an isotropic, homo- 
geneous, Gaussian random velocity field with zero mean, 
and have a short (ideally infinitesimal) correlation time 
r (Markovian approximation). Splitting B into a mean 



(large-scale) magnetic field (B) = B and a stochastic 
small-scale field b = B B, one derives the mean-field 
dynamo equation 0, 
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v IS 
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£ - 77V x B), 
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Here £ = (vxb) » a n B-ftVxB, is the turbulent EMF, 
where cto — — (t/3)(vt ■ V x vy) is the dynamo a-effect, 
proportional to the kinetic helicity and /3q — t(v^)/3 
is the turbulent magnetic diffusivity proportional to the 
specific kinetic energy of the turbulence. These equa- 
tions predict the exponential growth of the mean mag- 
netic field. One can also derive the equations for small- 
scale magnetic field correlations jl|, |2|, y| , which predict 
the exponential growth of small-scale fields on a shorter 
time scale. The kinematic theory then needs modification 
to take account of the back-reaction due to the growing 
Lorentz forces. 

In the quasi- linear approximation 0, , this is done 
by assuming that the Lorentz force induces an additional 
non-linear velocity component vjv, that is v = vo + + 
v^r, with vtv satisfying the perturbed Euler equation 

p{dv N /dt) = [B • Vb + b • VB]/(4tt) - Vp (3) 

and V • vjv = 0, where p is the fluid density and 
p the perturbed pressure including the magnetic field 
contribution. The turbulent EMF then becomes £ — 
(vt x b) + (vn x b), where the quasi-linear correction 
to the turbulent EMF £ N = (T(dv N /dt) x b). Here r 
is again a correlation time assumed to be small enough 
that the time-integration (over the correlation time) , can 
be replaced by simple multiplication. We ignore the (b 3 ) 
contributions to £, in the quasi- linear approximation, al- 
though these may indeed be negligible if the saturated 
small-scale field has a symmetric probability distribu- 
tion. One expects this approximation to give a reason- 
able estimate of non-linear effects, when the mean field is 
still weak, and be also analytically tractable. Some sup- 
port for the quasi-linear approximation also comes from 
EDQNM type closures of MHD 0]. 
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We will calculate £n in co-ordinate space represen- 
tation. We can eliminate the pressure term in vjv us- 
ing the incompressibility condition. Defining a vector 
F = a[B • Vb + b • VB], with a = r/(47rp), one then gets 



-N 



(F x b) - ([V(V~ 2 V-F)] x b), 



(4) 



where V~ 2 is the integral operator which is the inverse 
of the Laplacian, written in this way for ease of notation. 
We will write down this integral explicitely below, using 
— l/47rr to be the Green function of the Lapalacian. We 
see that £n has a local and non-local contributions. 

To calculate these, we assume the small-scale field to 
be statistically isotropic and homogeneous, with a two- 
point correlation function (6j(xi, i)6 J -(x2, t)) = Mij(r,t), 
where r — Xi — X2, r = |r| and 



Ma = M 



^-C-^)]+M L (^)+He ijf r f . (5) 



Mj,(r, t) and Mpj(r, t) are the longitudinal and transverse 
correlation functions for the magnetic field while H(r,i) 
represents the (current) helical part of the correlations. 
Since V • b = 0, M N = {l/2r)d{r 2 M L )/(dr). For later 
convenience, we also define the magnetic helicity correla- 
tion, N(r, t) which is given by H = —(N"+4:N'/r), where 
a prime / denotes derivative wirth respect to r. In terms 
ofb, wehaveM L (0,i) = (b 2 )/3, 2H(0,t) = (b-Vxb}/3 
and 2N{0, t) = (a • b) /3 (where b=Vxa). 

The local contribution to £n is easily evaluated, 



(F x b) = -aM L (0, t)(V x B) + 2aH (0, i)B (6) 



At this stage (before adding the non-local contribution) 
there is indeed a non-linear addition to the diffusion of 
the mean field (the — aMi(0, t)(V x B) term). Let us now 
evaluate the non-local contribution. After some algebraic 
simplification, this is explicitely given by the integral 



= 2e 



-«[V(V- 2 V-F)]xb)), 

d 3 rr 3 dM mk {v,t) dBi{y,t) 



ijk 



An r 3 



dr l 



(•7) 



where y = r + x. Note that the mean field B will in 
general vary on scales R much larger than the correla- 
tion length I of the small-scale field. We can then use the 
two-scale approach to simplify the integral in Eq. J7J). 
Specifically, assuming that (l/R) < 1, or that the varia- 
tion of the mean field derivative in Eq. J7J , over I is small, 
we expand dBi(y,t)/dy m , in powers of r, about r = 0, 



dBi 
dy m 



dBi 

dx m 



i 92 Bl 
lp dx m dxP 



d 3 B, 



dx m dxPdxi 



where we have defined n, = rj jr (we will soon see why we 
have kept terms beyond the first term in the expansion) . 



Simplifying the derivative dM m k(r,t)/dr l using Eq. (JSJ 
and noting that eijkfjrk = 0, we get 



dM m k 

dr l 



rjeijkl 



. (Ml — M J 



N 



n m hi + M' N niS 

ink 



+ He mk i +rH'n f nie mk f]. 



(9) 



We substitute © and JjJJ into Q, use j(dil/4Tr)n i n J = 
Sij/3, and J(dCl/in)ninjnkni = [5i 3 8 k i + 5 ik S 3 i + 
5uSj k ]/15 to do the angular integrals in Q, to get 



cNL 

- — 



-aM L (0,t)(V x B) + — N(0,t)V 2 B 
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drrM L (r, t) 



V 2 (V x B) (10) 



The net non-linear contribution to the turbulent EMF is 
£ N = £^ + £^ L 1 got by adding Eq. © and Eq. (TO). We 
see firstly that the non-linear diffusion term proportinal 
to V x B has the same magnitude but opposite signs 
in the local (Eq. ©) and non-local (Eq. (JTUJ)) EMF's 
and so exactly cancels in the net £n- This is the often 
quoted result 0, 0] that the turbulent diffusion is not 
renormalised by non-linear additions, in the quasi-linear 
approximation. However this does not mean that there is 
no non-linear correction to the diffusion of the mean field. 
Whenever the first term in an expansion is exactly zero 
it is neccessary to go to higher order terms. This is what 
we have done and one finds that £n has an additional 
hyperdiffusion £ HD — 7/hdV 2 (V x B), where 



2a f°° 

tjhd = -=- / dr rM L {r,t). 
5 Jo 



(11) 



Taking the curl of £n, the non- linear addition to the 
mean- field dynamo equation then becomes, 

V x & = K/ + /i M V 2 ]V x B - r/ ffD V 4 B (12) 

Here olm = a(h -V x b)/3 is the standard non-linear 
correction to the alpha-effect and Um — a(a. ■ b)/5 

is an additional higher order non-linear helical correction 
derived here. 

One can check that the hyperdiffusion coefficient t]hd 
is positive definite, by writing (jl 1|> in Fourier space. 
The longitudinal magnetic correlation is given in Fourier 
space by Mi(r,t) = 2/ dkEM(k,t)(ji(kr)/kr), where 
EM(k,t) is the magnetic power spectrum, and ji(x) the 
spherical Bessel function. Using this relation one gets 
Vhd = (4a/15) J dk(E m (k,t)/k 2 ), which is clearly posi- 
tive definite. The magnitude of t/hd ~ (2a/5)Mi(0, t)l 2 . 
So the importance of hyperdiffusion, relative to the 
turbulent diffusion is given by (i]hd/R 4 )/((3o/R 2 ) ~ 
(2/5)(b 2 /47r ( ov 2 )(^/i?) 2 . So for equipartition small scale 
fields, hyperdiffusion is only important in mean-field evo- 
lution, for moderate scale speparations l/R < 1. It could 
play an important role for example, in the "self-cleaning" 
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evolution seen in simulations of Brandenburg , by caus- 
ing a non-linear cascade of power from large-scale fields to 
nearby (in scale) smaller and smaller-scale fields. In case 
l/R <C! 1, the usual alpha-suppression gg, arising from he- 
licity conservation (and consequent growth of o>m of the 
right sign to cancel ao), is expected to lead to mean- field 
dynamo saturation, rather than hyperdiffusion. In both 
situations, since l/R is smaller than unity, non-linear cor- 
rections of higher order than hyperdiffusion are expected 
to be smaller by further factors of (l/R)- Analogous ef- 
fects are expected to be crucial for small-scale dynamo 
saturation, to which we now turn. 

The non-helical small-scale dynamo: It is well known 
that small-scale magnetic fields can grow under the ac- 
tion of a stochastic velocity field vy, even if the flow is 
non-helical, provided R m is greater than a critical value 
of order 100. The kinematic problem is well studied in 
the literature . We consider now how such a non- 
helical, small-scale dynamo saturates. For this one can 
neglect the subdominat effect of vo and also the B cou- 
pling to b (since b is expected to grow much faster than 
B). To model the effects of non-linearity, and in anal- 
ogy to the above quasi-linear treatment, we assume that 
Lorentz forces due to the growing small-scale field in- 
duces as additional non-linear velocity component vat, 
satisfying an equation analogous to Eq. ©. That is 
now v = + vat, with vjv = a[b ■ Vb — Vp], and 
V • vat = 0. Here once again p includes the perturbed 
magnetic pressure. Using the incompressibility condi- 
tion, one can again write v^r = v L + v NL , as the sum 
of a "local" term v L = ah ■ Vb = f, and the non-local 
"pressure" term v NL = -V(V~ 2 V • f). 

The stochastic Eq. can now be converted into the 
evolution equation for Ml- The detailed derivation of 
this equation with a different non-linear velocity compo- 
nent (modelled as an ambipolar type drift) is given in [3j. 
The major difference here is the form of the non-linear 
term, which we evaluate explicitely below. We get 



dt 



2 d 4 dM L , 



dr 



dr 



GM L + K 



(13) 



deal with fourth moments of the fluctuating field b. As 
in 0, we use a Gaussian closure to write these fourth 
moments in terms of the second order moments. The 
local contribution to K \j , involving only the v L terms in 
Eq. I|14|). can then be simply evaluated to give 



K^=2aM L (Q,t)V 2 M ir 



(15) 



The non-local contribution K^ L , involving the pressure 
term, v NL , is again expressible as an integral, after using 
the Green function of V 2 . The x and y-derivative terms 
in Eq. l|14f> give equal contributions and we get, 



k nl _ R (y) 
y jpq 



la d 3 u u p dM mq (u, t) dMu (X, t) ^ 



47T U 3 



du l 



dX' 



where X = u + y — x = u + R = u r. 

Since = K^+K^ L , one gets an integro-differential 
equation for the evolution of Ml, which is not analyti- 
cally tractable. One can however make analytic headway 
in two limits r = |x — y| 3> I, and r < 1, where l(t) is 
now the length scale over which Ml(t, t) is peaked. (For 

1/2 

example, during kinematic evolution, I = r^j ~ L/Rm , 
where L is the velocity correlation length 0, 0). F° r 
r > I, the integral (|16|l can then be evaluated by tak- 
ing the limit r = |x — y| u, and again expanding 
dMu(K, t)/dX' m , in powers of u, about u = 0. So 



0M K (X) dM H (R) 



d 2 M, 



dX J 



dR m 
u 2 n s nt 



dR r < 
d 3 M H 



dR s 



dR m dR s dR t 



(17) 



where we have defined now n, = Ui/u. Substituting the 
above expansion into Eq. (|16fl . using again Eq. J5j to 
simplify the derivative dM mq (u,t)/du l , and evaluating 
the angular integrals with the help of various moments 
of n, defined above, we get 



K% L (r,t) = -2aM L {Q,t)X7 2 M l3 



• 2t 7H £>V 4 Mtf. (18) 



where we have defined k = X] + ?l(0) — Tl(t), G — 
-2(T' L ' + AT'Jr) and T L (r) is the longitudinal correla- 
tion function of tvt defined analogous to Ml (see 0] ) . 
The diffusion k includes microscopic diffusion (77) , a scale- 
dependent turbulent diffusion (7l(0) — Tl(t)) and the 
G(r) term allows for the rapid dynamo generation of 
magnetic fluctuations by velocity shear 0, 0, . 

The non-linear contribution is K — (r^j /r )Kij where 



We see that, for r 3> /, again has a diffusion 

term which exactly cancels the corresponding term in 
Kh, leaving behind pure non- linear hyperdiffusion, that 
is Kij = Kt<+K$ L = -2n H D V*Af tf (r,t). (Also no odd 
derivative terms appear in the absence of hclicity). So 



Kir t) = -^R — 

r 4 Qr 



dr V r 4 dr dr 



(19) 



Ki 



^S(([<(y) + <%)Mx)My)>) 

+#2«bp (x) + ^(^(x^y))). (14) 



Here the operator R^ q = eu m e mp q(d / dx l ) and R^ q 



£iim£mpq{d / dy ) . In order to evaluate K, we need to 



Now consider the other limit r C 1. In this limit, 
one can assume r C « in Eq. I|16fl . and expand around 
r = 0. In fact since the first term (with r = 0) nei- 
ther vanishes, nor cancels exactly (see below), we 
need to only keep this term. Putting r = in Eq. i|16|) 
for Kf: L , straightforward but tedious algebra gives, 
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K^ L = 8Sij J(du/u) (M' L ) 2 . Again adding the local and 
non-local terms we get 

K(r,t) = 2aM L (0,t)^(r^) + 8 [°° ^(M' L f 

(20) 

Note that now for r <C I, K is still in the form of 
a non-linear diffusion, albeit with a partial cancella- 
tion due to the addition of the positive definite K NL 
contribution. One can check for specific forms of Ml, 
that the above K(r, t) does indeed lead to non-linear 
dissipation. For example, for a model M^(r, t) — 
Ml(0, t) exp(— r 2 jl 2 ), strongly peaked about I, we get 
K(0,t) = -24aMf(0,t)/Z 2 and even for r ~ 1/2, K is 
still negative ~ —6aM 2 (0,t) /I 2 . Of course as one goes 
to larger r ~ I, one has to keep higher order terms in the 
expansion around r = 0. 

From Eq. (|19|l and (|2U|) we see that the back reac- 
tion on the non-helical small-scale dynamo, due to the 
growing Lorentz force, can be characterised in this model 
problem, as non-linear diffusion for small r <C I (yet 
partially compensated by a constant), transiting to non- 
linear hyperdiffusion for r 3> I. And the damping of 
Ml in both regimes have damping coefficients which are 
themselves proportional to the magnetic energy density, 
or Mi(0,i). This means that as the small-scale field 
grows and Mi(0,i) increases, the non-linear damping 
would increase leading eventually to a saturated state. 
The properties of the saturated state requires detailed nu- 
merical solution, which we hope to return to elsewhere. 
But taking a clue from our earlier work where one 
had purely additional non-linear diffusion, the station- 
ary state could have a correlation function, which corre- 
sponds to a " ropy" small scale field, and with peak mag- 
netic fields of order the equipartion value. The way this 
is altered due to hyperdiffusion at larger r, albeit where 
Ml is subdominant, will be interesting to examine. 

Conclusion: We have examined here consequences of 
one popular model of non-linear back reaction on the dy- 
namo. For the mean-field dynamo, it has been thought 
that turbulent diffusion is not renormalised at all by the 
Lorentz forces. We have clarified that this is valid only 
at the lowest order, and at a higher order (to which 
one must go if the lower order term is exactly zero), 
one gets additional non-linear hyperdiffusion of the mean 
field [8j- Such hyperdiffusion may not be crucial for the 
mean field dynamo saturation, for I / R <C 1 . But it could 
have interesting consequences for how the field eventually 
"self-cleans" (or orders) itself during saturation. Further, 
when a similar model is applied to discuss the saturation 
of the non-helical small-scale dynamo, one obtains an 
intriguing combination of non-linear diffusion and hyper- 



diffusion, which governs how the small-scale fields reach a 
saturated state, due to the back-reaction of the Lorentz 
force, ft remains to solve the above equations numeri- 
cally and also elucidate the conditions when other possi- 
ble models |3|, for dynamo saturation are applicable. 

I thank Axel Brandenburg and Anvar Shukurov for 
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